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Complexity of modern power network and Large disturbance results voltage 
collapse. So, voltage security analysis is important in power system. 
Indicators are helpful in voltage stability analysis, as they give information 
about the state of the system. In this paper a new indicator namely 
Distribution System Stability Indicator (DSSI) has been formulated using the 
information of Phasor Measurement Unit (PMU).The proposed indicator 
(DSSI) is tested on standard IEEE 33 bus radial distribution system. The 
suggested indicator is also applicable to the equivalent two bus system of a 
multi-bus power system. The proposed indicator is calculated for different 
contingent conditions at different system load configurations. The result of 
DSSI is verified with the standard indicator (VSI) which proves applicability 
of the proposed indicator. The bus voltages of all the buses at base loading 
and at maximum loading are evaluated for base data and for tripping of most 


critical line. 
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1. INTRODUCTION 

Radial Distribution system has a large number of buses. So, it is not possible to maintain adequate 
bus voltage at each bus. This will increase the chance of voltage instability at the buses [1]. Increase in load 
demand will reduce the voltage at the receiving end which increases reactive power demand, if the power 
system is not able to supply the same there will be voltage collapse [2]. For proper operation of power system 
the voltage stability of the power system should be increased [3]. To improve power system stability and 
security, planning and operation of power system is essential [4]. Stability analysis may be of two forms- 
static and dynamic. In static analysis load flow based simulation is performed and in dynamic analysis time 
domain based analysis is performed [5]. With the implementation of phasor measurement unit (PMU) in 
power system for voltage security analysis the security level of power system has been increased. Online 
based methods and offline based methods both have been studied by the researchers for gaining information 
about voltage stability of power system using phasor measurement unit [6]. 

Different techniques have been used for voltage stability analysis. P-V curve and Q-V curve 
method [7-9] are the traditional methods of voltage stability analysis. Many indices such as SVSI [10], 
L-index [11], Lin and VCPI [12], integrated voltage stability indicators [13] etc. can predict voltage 
instability in a power system. Other methods such as radial basis function and kohonen’s self-organising 
feature map can also be found in literature for voltage stability analysis [14], [15]. In recent years a lot of 
work has been done in the area of voltage stability. Voltage stability analysis was done by the researchers 
using nodal analysis [16], modified differential evolution [17], coupled single port circuit [18], circuit theory 
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approach [19], catastrophe theory [20], chaotic particle swarm optimization algorithm [21], vector regression 
model [22] etc. Along with the other methods phasor measurement unit based schemes are also been used in 
power system to monitor voltage stability [23], [24]. PMU in substation gives voltage and current phasors 
information in a microsecond when the reading is taken. It capable the system for online monitoring. By 
using PMU data in the current literature accurate rating can be obtained. 

This paper deals with a new voltage stability indicator which is used for monitoring voltage stability 
condition in base case as well as at maximum loading condition considering different line tripping condition. 
The developed distribution system stability indicator (DSSI) is tested on phasor measurement unit 
incorporated power system model of standard IEEE 33 bus radial distribution system. The result of DSSI 
indicator can accurately identify the critical conditions of line outage. The result is verified with the standard 
VSI indicator [25].The change of bus voltages from base value to tripping of Line Number 23 is also shown 
for base load and maximum load. 


2. RESEARCH METHOD 

In this paper phasor measurement unit (PMU) incorporated voltage security analysis has been 
performed. PMU is allocated in every bus to give both magnitude and angle information of bus voltage and 
current. 

Figure 1 shows the flowchart of the proposed method. 


Voltage and Current 
measurement by 


PMU 


Calculation DSSI of the 
system by tripping 
different lines 


Identification of Critical 
Line Tripping cases 


Voltage Security 
Analysis of the system 


Figure 1. Flowchart of the proposed method 


By using PMU at every bus of the concerned multi-bus system the voltage and current can be 
measured accurately at the same instant of time. These real time data from PMU are used in the proposed 
method to calculate the distribution system stability indicator (DSSI). Then this technique find the value of 
DSSI indicator for different conditions and from them identify the critical cases. From the results the voltage 
security analysis is performed. 


3. DERIVATION OF THE PROPOSED INDEX FROM GAUSS-SEIDEL EQUATION 

Figure 2 depicts a simplified two bus network where source is connected with bus number | (slack 
bus) and a load is connected to bus number 2 (load bus). To make the system observable in respect to 
voltage, phasor measurement unit is installed in both the buses. The current flow between bus number 1 and 
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bus number 2 is I. The load absorbs apparent power S2 = P2 + j Q» from the system. Two buses are connected 


together by impedance” Zag T Reg + jE eq À 
Assuming slack bus voltage equal to unity, 


_ Peen — Proad z Qcen — Qiroaa 


Reg = 
q 2 2 1“ eq 2 2 
Pen + QGen Foen + Qen 


V,20 V,26 


Figure 2. PMU incorporated simplified two bus network 


The voltage at the receiving end bus that is at Bus number 2 is given as 


Vp = V, — acon een) (1) 


The Gauss-Seidel load flow Equation for the two bus system can be written as 
— 1 [P2-jQ2 _ 
n = EP a] (2) 


Where, V, is the voltage at bus number 1, Vz is the voltage at Bus number 2 and Y» is the 
admittance of Bus number 2. 
Equation (2) can be written as 


oo 
-a YV, = VÝ = h (3) 


But, 1, =1, =1 (4) 
From the two bus system of Figure 2 current I, can be written as, 


l= net (5) 
eq 


Equating Equation (3) and Equation (5) it can be written that, 


V1Z0-V2265 _ P2-jQ2 
= * a Yi2V1 
Zeq v2 
V4 Z0—-V226 P2-jQ 
=+ Y2 V1 = -= Z z 
Zeq v2 
V140-V24ő P2-j 
= Bee YV; = 2—JQ2 
Zeq V22-6 
= V{20—-V225+¥12ViZeq _ P2-jQ2 
Zeq E V22-6 
2 _ : 
=> Vi V24 — ô — V3 + Zeq¥12ViV24 = ô = ZeqP2 —J Q2Zeq (6) 


Taking real part from Equation (6) it can be written as 
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VaV, cos 8 — V? + Yi2VjV2Zeq COS Ô = ZeqP2 
= ZeqP2 — Yı2VıV2Zeq cos 5 = VıV; cos ô — V? 


_  WWVzcos 8-v2 
eq ~~ Pa- Y12V1V2 cos 8 


LZ (7) 


Taking derivative of both sides of Equation (7) with respect to V, it can be written that, 


d Zeq _ (Vı cos 5—2V2)(P2—Y12V1V2c0s ô)—(V1V2 cos ô— v2)(0-Y12 Vıcos ô) 
dV2 (P2— ¥12V1V2 cos 8)? 


(8) 


After simplification Equation (8) can be written as, 


dZeq _ Y42V1Vf cos 6—2P2V2+ P2V: cos ô 
dv, (P2— Y42V4V2cos 5)? 


(9) 


At critical point of voltage stability, 


d Zeq 
dv, 


So, Equation (9) it can be written as, 
(Y%42V,cos 5)V2 = (2P2)V2 + P Vi COS ô — 0 
So, critical voltage at the receiving end of the equivalent two bus network is written as, 


_(p2_ 3 2g) 05 
Vane [P2-(P}-Y12Vf P2cos 75) `] (10) 


Y12V,cos 6 


Distribution system stability indicator (DSSI) can be formulated by the difference between the 
receiving end voltage of the equivalent two bus distribution system at base case and critical voltage at the 
receiving end of the equivalent two bus distribution system and can be written as, 


DSSI = Vo — Voeri 


Where, V> is the receiving end voltage of the two bus equivalent system and V>,,; is the critical 
voltage at the receiving end of the same system. So, at the critical point of voltage stability the difference 
between these two will be zero. The value of the proposed indicator will provide information about the 
voltage stability of the system. The higher value of DSSI indicator means the system is less secure in terms of 
voltage stability and lower value of DSSI indicator means more secure is the system in terms of voltage 
stability. 


4. SIMULATION AND RESULTS 

The program is simulated in MATLAB software and the developed indicator (DSSD is implemented 
in standard IEEE 33 bus radial distribution system. The system is consist of active power load of 3715 kW 
and reactive power load of 2300 kVA. The whole system is reduced into two bus equivalent network. The 
result of DSSI indicator is compared with standard VSI indicator. The simulation results are described in the 
following sections. 


4.1. Identification and Verification of Critical Contigency Cases 

Contingencies are the harmful situations on power system that may occur during steady state 
condition of power system. The contingency condition may be outage of generators or tripping the line of a 
power system. To check the accuracy of the proposed indicator, DSSI is calculated for different contingent 
conditions of the said IEEE 33 bus radial distribution system. The results are verified by the standard VSI 
indicator. 

Table | shows the results of the test. 
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Table 1. Contingency test results at base loading condition 


System System DSSI VSI 
Configuration Configuration 

number 
1 Base case 0.9397 0.7332 
2 Line 17 tripped 0.9398 0.7333 
3 Line 19 tripped 0.9399 0.7342 
4 Line 24 tripped 0.9401 0.7348 
5 Line 30 tripped 0.9403 0.7357 
6 Line 29 tripped 0.9409 0.7381 
7 Line 23 tripped 0.9412 0.7395 


Figure 3 shows the graph of DSSI indicator and VSI indicator for base case and various contingency 
conditions at base loading depicted in Table 1. 
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Figure 3. Plot of (a) DSSI and (b) VSI indicator for different configurations at base loading 


Table 1 and Figure 3 clearly show that most severe contingency is system configuration 7 i.e. 
tripping of line number 23 because the value of the DSSI indicator is highest in this case. Proposed indicator, 
DSSI gives accurate result which is verified by the VSI indicator. 

After performing contingency test at base load the loads of all the buses have increased successively 
with an increment of 5% until the load flow program fails to give solution. Maximum loading at which load 
flow program runs is known as the critical loading of the system. At maximum loading condition 
contingency test is performed to evaluate the change in the DSSI indicator value of the concerned power 
system. 

Table 2 shows the result of contingency test at maximum loading condition. 


Table 2. Contingency test results at maximum loading condition 


System System Configuration DSSI VSI 
Configuration 

number 
1 Base data at maximum loading 0.9755 0.8877 
2 Line 17 tripped 0.9756 0.8878 
3 Line 19 tripped 0.9757 0.8881 
4 Line 24 tripped 0.9758 0.8884 
5 Line 30 tripped 0.9759 0.8889 
6 Line 29 tripped 0.9760 0.8899 
7 Line 23 tripped 0.9761 0.8906 


Figure 4 shows the graph of DSSI indicator and VSI indicator for base case and various contingency 
conditions at maximum loading point. 
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Figure 4. Plot of (a) DSSI and (b) VSI indicator for different configurations at maximum loading 


Table 2 and Figure 4 proves that at maximum loading condition the worst situation is system 
configuration number 7 i.e. tripping of line number 23. The accuracy of the DSSI indicator is checked with 
VSI indicator. 

At base loading, the effects of tripping of various lines affect the state of the system and accordingly 
the value of the indicator changes from base case. So, the change in the indicator values will give indication 
about tripping of line. With the change in load demand, the indicator value changes at every contingency. At 
the maximum loading condition the value of the indicator will be different from the indicator values at base 
loading for the same contingencies. 

Figure 5 displays the change of DSSI indicator from base loading to maximum for different 
contingencies. 


CONFIGURATION NUMBER 
—+— DSSI at base loading —E— DSSI at maximum loading 


Figure 5. Comparison of DSSI value at base loading and at maximum loading 


Figure 5 shows that at maximum loading condition the value will be higher than those obtained at 
base loading condition. 


4.2. Effect of Contingency on Bus Voltages 

Voltage magnitudes at different buses is also changes due to contingency. Figure 6 compares the bus 
voltages for base case and for tripping of Line number 23 at base loading. 

Figure 6 shows that for base data and for tripping of line number 23 there is a small change in bus 
voltages. 
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At the point of maximum loading the system is more prone to voltage instability and the bus 
voltages may be less than the allowable limit. So, it is essential to check the bus voltages at this stage of 
power system. With the increase in loading the bus voltages will decrease. 

Figure 7 shows the bus voltages at maximum loading for all the buses of the concerned distribution 
system for base data and for outage of line number 23. 


—+—Base data —w— Line 23 out 


Figure 6. Comparison of original bus voltages and bus voltages for outage of line number 23 at base loading 


=e Base ditam Line 23 out 


Figure 7. Comparison of actual bus voltages and bus voltages for tripping of line 23 at maximum loading 


Figure 7 shows that the variation in bus voltages with actual data and for tripping of Line number 23 
at maximum loading is very small. 


5. CONCLUSION 

By applying the proposed indicator the critical contingency cases of the system can be identified 
very easily which is also verified by the standard VSI indicator. At maximum loading the security state can 
be accessed by the values of DSSI indicator for various line tripping conditions. If the bus voltages are below 
the adequate value set for the system, control action can be taken to improve the bus voltages and voltage 
stability for the system can be improved. From the knowledge of the indicator the power system operator will 
be able to check the voltage stability limit. The DSSI value will also able the user to predict the contingency 
cases. So, it will be possible to provide protection to the power system. 
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